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ABSTRACT 

HST V and I-band observations show that the gravitational lens B1359+154 consists 
of six images of a single Zg = 3.235 radio source and its star- forming host galaxy, 
produced by a compact group of galaxies at ~ 1. VLB A observations at 1.7 GHz 
strongly support this conclusion, showing six compact cores with similar low-frequency 
radio spectra. B1359-I-154 is the first example of galaxy-scale gravitational lensing 
in which more than four images are observed of the same background source. The 
configuration is due to the unique lensing mass distribution: three primary lens 
galaxies lying on the vertices of a triangle separated by 0'.'7 ~ 4/i~^ kpc, inside the 
l'!7 diameter Einstein ring defined by the radio images. The gravitational potential 
has additional extrema within this triangle, creating a pair of central images that 
supplement the "standard" four-image geometry of the outer components. Simple mass 
models consisting of three lens galaxies constrained by HST and VLBA astrometry 
naturally reproduce the observed image positions but must be finely-tuned to fit the 
flux densities. 
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1. Introduction 

The vast majority of the ~ 60 known arcsecond-scale gravitational lens systems consist of two 
or four detectable images, consistent with the generic lensing properties of smooth, isolated and 
centrally steep mass distributions (e.g. Blandford & Kochanek 1987). There are very few cases in 
which a "non-standard" number of images may have been observed.^ APM08279-I-5255 (Ibata et 
al. 1999) and MG1131-I-0456 (Chen & Hewitt 1993) each contain a central component that could 
be an additional lensed image created by a sufficiently large galaxy core or shallow mass profile 
(Narasimha, Subramanian & Chitre 1986; Rusin k Ma 2000). Alternatively, APM08279+5255 
may be a special class of imaging produced by an edge-on disk (Keeton &: Kochanek 1998), 
while the central component in MG1131+0456 could be weak AGN emission associated with the 
lensing galaxy. In addition, MG2016+112 (Lawrence et al. 1984; Garrett et al. 1996) exhibits a 
rather complicated image morphology consisting of four primary components, one of which has 
three subcomponents in VLBI maps. This may be a compound lens produced by two galaxies at 
different redshifts (Nair & Garrett 1997). Benitez et al. (1999) suggest an alternative model, but 
it might not fully explain the VLBI observations. 

Complex interplay between mass distributions can lead to lens systems with more than four 
images of a source. For example, ellipsoidal mass distributions perturbed by shear fields may 
produce configurations with six or eight images arranged about the tangential critical curve 
(Keeton, Mao & Witt 2000b). However, this requires that the relative magnitudes and orientations 
of the internal and external shear axes be finely tuned, and the resulting cross-sections are quite 
small. Compound mass distributions have been shown to be far more efficient at producing 
a variety of complex geometries in which five or more images may be formed over a range of 
radii (Kochanek &; Apostolakis 1988). This has been observed in the cluster-lensing case, where 
substructure in the gravitational potential created by the resident galaxies can qualitatively alter 
the lensing properties that one would expect for a smooth halo mass distribution (e.g. Natarajan 
et al. 1998; Meneghetti et al. 2000). A dramatic example of this is CL0024+1654, which exhibits 
eight images of a single blue background galaxy (Kassiola, Kovner & Fort 1992; Wallington, 
Kochanek & Koo 1995; Colley, Tyson & Turner 1996; Tyson, Kochanski & Dell' Antonio 1998). 

Early- type galaxies preferentially participate in lensing (Kochanek et al. 2000a), so a 
large fraction of lens galaxies should be members of small groups and clusters due to the 
morphology-density relation (Dressier 1980; see also Keeton, Christlein &: Zabludoff 2000a). 
Indeed, many lens systems are known to be perturbed by nearby galaxies (e.g. B2319+051; Rusin 
et al. 2000b), groups (e.g. PG1115+080; Schechter et al. 1997), or clusters (e.g. QSO 0957+561, 
Fischer et al. 1997; RXJ0921+4529, Muiioz et al. 2000). A handful of systems are lensed by 
more than one primary galaxy (e.g. B1127+385, Koopmans et al. 1999; B1608+656, Koopmans 



*B1933+503 (Sykes et al. 1998) and B1938+666 (King et al. 1997) contain ten and six lensed radio components, 
respectively, but this is due to the imaging of a multi-component source rather than any exotic properties of the 
lensing potential. 
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&; Fassnacht 1999), but in these cases standard image geometries are produced despite merged 
caustics. Several additional lenses are observed to have faint satellite galaxies as companions 
inside or near the Einstein radius (e.g. MG0414+0534, Schechter & Moore 1993; B1030+074, 
Xanthopoulos et al. 1998; Bl 152+199, Rusin et al. in preparation). While the likelihood of finding 
nearly equal mass galaxies close enough to have merged caustics is predicted to be only ~ 1% 
(Kochanek &; Apostolakis 1988), the probability of finding fainter satellites near the primary lens 
is not small because galaxy luminosity functions diverge for faint galaxies - all lenses should have 
faint neighbors, as discussed in the Appendix. While such systems typically produce regions of the 
source plane in which more than four images can form, the small size of the companion galaxies 
means that the cross-sections for creating non-standard image geometries are not significant. 

The gravitational lens system B1359+154 (Myers et al. 1999), discovered in the Cosmic Lens 
All-Sky Survey (CLASS; e.g. Myers et al. 1995), has been suspected of containing more than four 
lensed images. Observations with the Very Large Array (VLA) and Multi-Element Radio-Linked 
Interferometer Network (MERLIN) show a total of six radio components (Myers et al. 1999; Rusin 
et al. 2000a). Four of these components (A-D) are arranged in a typical quad-lens configuration 
(maximum image separation of 1'.'7), with two additional components (E and F) residing within 
the ring defined by the outer images. Preliminary radio spectral studies of B1359+154 at high 
frequency (Myers et al. 1999) suggested that E had a slightly flatter spectrum than A-D, and 
therefore that the central components may be core-jet emission associated with a weak AGN in 
the lensing galaxy or galaxies, as in the case of B2045-I-265 (Fassnacht et al. 1999). When the 
spectra were extended down to 5 GHz, however, there appeared to be less disparity among the 
radio components (Rusin et al. 2000a). Subsequent VLA observations at 15 GHz have failed 
to decisively confirm that E is fiatter than A-D at high-frequency, or detect component F. 
Spectroscopy with the Keck II telescope determined the source redshift to be Zg = 3.235 (Myers et 
al. 1999). Adaptive optics observations of B1359-I-154 conducted with the Canada-France Hawaii 
Telescope (CFHT) in the infrared K' band (Rusin et al. 2000a) detected counterparts to radio 
components A-D, and discovered three extended emission peaks (K1-K3) bracketing the expected 
positions of E and F. K1-K3 were identified as three possible lensing galaxies, comprising the core 
of a compact galaxy group. Evidence of an arc connecting A, B and C was also observed, along 
with a weaker emission feature associated with component E. The compound defiector system not 
only explained why attempts to model the outer four components using a single galaxy had failed 
(Myers et al. 1999), but offered the means of creating additional extrema in the lensing potential. 
This opened the possibility that at least one of the central components is a lensed image. 

In this paper we present powerful new evidence from observations with the Hubble Space 
Telescope (HST) and Very Long Baseline Array (VLBA) that B1359-I-154 consists of six images 
of a single background source, lensed by a compact group of galaxies at ~ 1. In §2 we 
present VLBA observations of B1359+154 and investigate the low-frequency radio spectra of the 
components. In §3 we discuss and analyze HST V and I-band observations, which offer compelling 
evidence for the six-image hypothesis. In §4 we use preliminary mass modeling to demonstrate that 
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B1359+154 can be naturally explained as a true six-image lens system. Finally, in §5 we discuss 

the prospects for obtaining improved constraints on the lensing mass distribution of B1359+154, 
and ultimately using the system to study the structure of small galaxy groups at high redshift. 

2. VLBA Imaging of B1359+154 

VLBA observations of B1359+154 were obtained at 1.7 GHz on 1999 December 10 for an 
on-source integration time of 2 hr, and again on 2000 Aug 28 for 8 hr. These observations 
were performed by iterating between the target source for 3 min and the nearby phase-reference 
calibrator B1413-I-150 for 1.5 min. The data for each of the two epochs were calibrated separately 
within the AIPS data reduction package, then mapped and model-fit using DIFMAP. The flux 
densities of the model-fit components exhibited little variation 3%) between epochs, so the two 
data sets were combined to improve the signal-to-noise. The resulting naturally-weighted map 
of B1359-I-154 is presented in Fig. 1, and has an rms noise level of 45 //Jy/beam. All six radio 
components seen in previous VLA and MERLIN observations are easily detected by the VLBA, 
and maps of the individual components are shown in Fig. 2. Components A, B, and C each 
exhibit a compact core with associated jet emission. Lensing-induced parity reversal is evident in 
the relative orientations of the subcomponents. Components D, E and F arc each unresolved. The 
data were modeled within DIFMAP using a total of 9 gaussian components (Table 1): 6 compact 
cores (Al-Fl), 1 extended jet (A2) and 2 compact subcomponents (B2, C2). 

The compact nature of radio components E and F strongly argue for their identification as 
additional lensed images, rather than weak core-jet emission associated with the lensing mass 
(Myers et al. 1999). First, it would be rare for jet emission to remain unresolved at ~ 10 mas 
resolution. All but one pair of radio components shown to be similarly compact in CLASS VLBA 
follow-up observations is a pair of lensed images, the only exception being the binary quasar 
B0827-1-525 (Koopmans et al. 2000). There is also no evidence in the VLBA data to suggest that 
any emission bridge might be connecting E and F. Second, it is unlikely that both E and F could 
mark the cores of independent AGN within the deflector. One would expect the positions of AGN 
to be correlated with brightness peaks of the lensing galaxies, and this is not the case in the CFHT 
data. Third, the structure of all six radio components are morphologically consistent with a single 
background source. Specifically, if the source consists of a core and weak jet, it is likely that 
extended emission would be visible in the more magnified lensed images while the fainter images 
show only the compact core. This is exactly what is observed in B1359+154: A-C are bright and 
exhibit weak subcomponents; D- F arc faint and share identical unresolved morphologies. 

Next we construct low-frequency radio spectra of the components by combining the 1.7 GHz 
VLBA flux densities with those previously obtained at 5 GHz using MERLIN (Rusin et al. 2000a) 
and at 8.5 GHz using the VLA (Myers et al. 1999). One need not be concerned about combining 
data from different epochs as there is little evidence for any significant variability in the lensed 
source. The resulting radio spectra are plotted in Fig. 3 and exhibit striking similarities. Each of 
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Fig. 1. — VLBA 1.7 GHz naturally-weighted map of B1359-M54. The data is restored with a beam 
of size 13.9 x 12.5 mas at PA = —3.34°. The map rms noise level is 45 //Jy/beam. 
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Fig. 2. — VLBA 1.7 GHz naturally-weighted maps of B1359-I-154 components. Top: A, B, C. 
Bottom: D, E, F. The beam is 9.88 x 5.54 mas at PA = 1.11°. Each box size is 250 x 250 mas. 
Contours are at 3, 6, 9...X the map rms. 
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Fig. 3.— Radio spectra of B1359+154 components. Data from the 1.7 GHz VLBA, 5 GHz MERLIN 
(Rusin ct al. 2000a) and 8.5 GHz VLA (Myers et al. 1999) observations. We assume flux density 
uncertainties of 5% to account for source variability and model-fitting errors. 
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the six components has a spectrum that is rather flat (af y ~ —0.3, where oc u°^) between 1.7 
and 5 GHz, then falls off^ sharply (ccg^g ~ —1.2) between 5 and 8.5 GHz. The 1.7 - 5 GHz spectral 
indices for all six components (Table 1) agree within measurement errors. Even if the sharpness 
of the 5 GHz spectral break were somehow exacerbated by variability, it would still require that 
all six components vary in the same manner - something that would be impossible unless E and 
F are related to A-D. Therefore, the morphologies and low-frequency spectra of the B1359+154 
radio components strongly suggest that they are six images of a single background source. 

3. HST Imaging of B1359+154 

Using the WFPC2 instrument on HST we obtained images of B1359+154 in the F555W 
("V", 5200 sec) and F814W ("I", 5000 sec) bands on 2000 July 10. Each observation was divided 
into four subcxposures, at two dither positions, which were used to remove cosmic rays and bad 
pixels. Fig. 4 shows the resulting images. Optical counterparts to all six radio components are 
observed, along with extended emission from the three galaxies. The images were modeled as 
a combination of point sources and de Vaucouleurs models following the modeling procedures 
of Lehar et al. (2000) for CASTLES observations of other lenses. Astrometric and photometric 
results are presented in Table 2. The positions of images A-F derived from the HST data match 
the radio positions to an accuracy of 0'.'02, consistent with our internal error estimates. 

All six optical counterparts of the radio components are blue, with an average color of 
V"I= 0.5 mag and a dispersion in the colors of only 0.1 mag. The flux ratios of the 6 radio 
components and their optical counterparts differ by less than a factor of two between 1.7 GHz and 
the V-band, so they must be six lensed images. Three of the counterparts (A, B and C) are clearly 
extended in both the V and I-band images, consistent with the presence of an Einstein ring in the 
Rusin et al. (2000a) infrared image. The rest frame wavelengths of the V and I-band fllters at the 
source redshift of Zg = 3.235 are 1300A and 1900A , respectively, so the host galaxy of the radio 
source must have a strong ultraviolet continuum to produce the arcs. While we lack images in 
bluer filters, and the Myers et al. (1999) spectrum of the source is contaminated by AGN activity, 
the host galaxy appears to be a star-forming Ly-break galaxy. The emission lines from the source 
are relatively narrow (Lya and NV are well resolved) and lines characteristic of both AGN (Civ) 
and star formation (Cm]) are seen in the spectrum. 

The optical counterparts (G, G', G") of the three infrared sources K1-K3 are extended and 
red. We have relabeled the sources using the standard CASTLES convention in which G is the 
most luminous lens and G'/ G" are of lower luminosity. All three are seen in the I-band image and 
essentially vanish in the V-band image, with a color of V-I^ 3 mag. They are clearly three distinct 
galaxies which arc well modeled by de Vaucouleurs profiles, rather than compact star- forming 
regions in a larger galaxy. With separations of only Ah~^ kpc (for ~ 1.0 in an = 0.3 flat 
cosmological model), they are likely to be the core of a compact group. Unlike other lenses with 
secondary lens galaxies, the three galaxies have similar luminosities. Redshift estimates based on 
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Fig. 4. — HST I-band (left) and V-band (right) images of B1359+154. The top panels show the 
final images with the components labeled as in Table 2. The bottom panels show the residuals 
after subtracting counterparts to the six radio cores and the three lens galaxies. Crosses mark the 
centers of the subtracted components. 
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their optical colors and the mass estimates derived from our models of the system (§4) and the 
fundamental plane (Kochanek et al. 2000a) are mutually consistent. We find Zi = 1.35 ± 0.16, 
0.88 lb 0.06 and 0.94 it 0.07 for galaxies G, G', and G", respectively. The average redshift is 
zi = 1.05 zb 0.20 and the weighted average is zi = 0.94 it 0.04, but for simplicity we will just assume 
zi = 1 for future estimates. The scatter of the individual estimates is consistent with the scatter 
found by Kochanek et al. (2000a) for lenses with known redshifts. Such a group of galaxies is 
unlikely to be a projection effect given the separations and the size of the overall lens sample. 
Based on the results in the Appendix we estimate that the a posteriori probability of finding such 
a chance projection in a sample of 60 lenses is ~ 10^^. Wc conclude that the lens is the core 
of a compact galaxy group (see Hickson 1997 for a review of compact groups). With three lens 
galaxies of comparable luminosity and redshift inside their joint Einstein ring, the gravitational 
potential should have additional extrema inside the triangle formed by the galaxies. While images 
A-D resemble a standard four-image lens morphology, two additional images (E and F) lie in the 
region between the three primary lenses. Therefore, HST observations offer powerful evidence that 
B1359+154 is the first example of a true six-image galaxy-scale gravitational lens system. 

We also cataloged the nearby galaxies in the HST images using the SExtractor package (Bertin 

6 Arnouts 1996). The software classified the optical objects, and computed total magnitudes 
using Kron-type automatic apertures (see Lehar et al. 2000). Colors were computed using fixed 
apertures scaled to the F814W size of each object. Table 3 lists all the galaxies detected by HST 
within 20" of the lens, after visually confirming the SExtractor classification. The galaxies were 
labeled G^,^ in order of decreasing I-band brightness, or for the galaxies found in the Rusin 
et al. (2000a) infrared observations which were not detected in the I-band image. Table 3 also 
includes an estimate of the shear perturbation each neighboring galaxy can produce on the lens 
assuming they lie at the same redshift and have the same mass-to-light ratios (see Lehar et al. 
2000). The sum of the shear contributions is 7 = 0.074 at 0^ = 83°, although the largest source 
of shear, galaxy Gl, is probably a foreground galaxy. The contribution of all the other galaxies is 

7 = 0.049 at 9^ = -60°. 



4. Preliminary Mass Models 

Constructing any believable model for the lensing mass of B1359-I-154 is a significant 
challenge. At a minimum the model must include mass distributions to represent each of the three 
lensing galaxies plus environmental contributions to the lensing potential, cither in the form of 
external shear or a separate group halo. Since even the simplest models require a large number 
of free parameters, the best hope for producing a realistic mass model may ultimately rely on 
constraints from deep HST imaging of the extended arc emission (Kochanek, Keeton & McLeod 
2000b). For now our primary objective is to demonstrate that a six-image lens system can be 
naturally produced by the observed lensing galaxies. 

All models are constrained with 12 coordinates and 5 flux density ratios of the lensed radio 
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components, as derived from the VLBA data listed in Table 1. Only the core components Al-Fl 
are used at this time. We set a fit tolerance of 1 mas for the image positions and assume a 10% 
uncertainty on the flux densities. The galaxy positions are constrained according to the astrometric 
error bars given in Table 2. The goodness-of-fit parameter is therefore Xtot = Xpos + Xjix + Xgai^ 
which includes contributions to the fit from the positions of the radio components, (Xpos)^ their 
flux iXfix)' the positions of the lens galaxies {xiai)- '^^^ redshift of the group is set at = 1, 
and a flat cosmology with = 0.3 is assumed for all calculations. The models are optimized 
using the lensmodel software package (Keeton 2000) .0 



4.1. A Demonstration of Principle 

Smooth gravitational potentials must produce an odd number of lensed images (Burke 1981). 
However, images can be trapped and demagnified in the singular or near-singular cores of lensing 
galaxies (Narasimha et al. 1986). We can label the possible image configurations by A/B, where 
A is the true number of images and B is the observed number of images. Since A is always 
odd, B1359-I-154 must be either a 7/6 or 9/6 system. A single lens can generate a 7/6 image 
configuration, but it should consist of 6 images near the Einstein radius of the system (see Keeton 
et al. 2000b). The central location of images E and F and the comparable luminosities of the 
three galaxies (1.0:0.4:0.4 for G:G':G") strongly suggest that all three lenses are important and 
that we are looking at a 9/6 lens system. 

We begin by considering a compound mass distribution consisting of three singular isothermal 
spheres (SIS) in an external shear field. This model is unlikely to fit the data very well as it 
ignores galaxy ellipticity, but we use it to demonstrate the principles of compound deflectors in 
the context of B1359-I-154. From the expected caustic structure alone one can argue that 9/6 
systems are possible. Each SIS deflector has an associated radial caustic on the source plane, at 
which an inward-crossing source produces a pair of images at the corresponding galaxy center: one 
observable and one trapped in the core. The asymmetry of the potential will also lead to at least 
one tangential caustic, at which an inward-crossing source produces a pair of observable images 
at the corresponding critical curve. Therefore, a 9/6 system is formed when a source sits within 
one tangential and three radial caustics. Because three positive parity images will be trapped 
in the galaxy cores, four of the observed images must have negative parity (saddle points in the 
time delay surface) to ensure that the parities of all nine images sum to 1 (Burke 1981; Blandford 
&L Narayan 1986). Images A-D are arranged in a typical quad configuration, so we expect their 
parities to alternate in the usual way: A(-|-), B(-), C(-|-), D(-). This means that both E and F 
must correspond to saddle points (-), and we model the system under this assumption. 

The critical curves, caustics and time delay surface for the best-fit 3SIS model are plotted in 



'The code is publicly available at 



http://penedes.as.arizona.edu~ckeeton/gravlenf 
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Fig. 5. Because the separation of the galaxies is small compared to their Einstein radii, a single 
tangential critical curve encloses all three galaxies (Fig. 5a). This maps to the highly distorted 
tangential (astroid) caustic on the source plane (Fig. 5b). The large circular features in Fig. 5b 
correspond to the radial caustics associated with each SIS. These not only significantly overlap 
each other, but also overlap much of the astroid caustic creating a substantial 9/6 region in 
which the source resides. The cross-section of the model is dominated by 3/2, 5/3, 7/5 and 9/6 
geometries, but it can also produce 7/4 and 11/8 configurations. The 11/8 region is due to the 
tangential caustic folding over on itself and creating a swallowtail catastrophe. The probability of 
observing 7/4, 7/5, 9/6 and 11/8 systems will be greatly enhanced through magnification bias. 
Despite the complexity of the caustic structure, the topology of the virtual time delay surface 
(Fig. 5c) is quite generic (Schneider 1985; Blandford & Narayan 1986). Images A-D have the 
standard topology for four- image lenses. Images A-C are associated with a lemniscate critical 
contour encompassing two minima (A and C) and the saddle point (B). Image D is a saddle point 
associated with a limagon contour which would usually encompass images A-C in the large loop 
and the core of the lens galaxy in the smaller loop. In this case, however, the smaller loop contains 
all three lens galaxies. Images E and F lie at the saddle points of two new lemniscate critical 
contours with maximum-saddle-maximum topologies, where the maxima sit on the galaxy cores. 

The 3SIS model has a best-fit of = 1570 for NDF = 10, and its parameters are listed in 
Table 4. The model poorly fits the positions of images A (Xpos.A — 323), B (Xpos,B — 659) and F 
{xl^^^p = 299). The flux density of image A is also significantly lower than observed (predicted 2.4 
™Jy5 Xjix A ~ ^'-*)- Given the simplicity of the model, the errors in the image positions are of little 
concern and should be greatly improved by adding ellipticity or adjusting the radial mass profiles 
of the galaxies. The mismatch of the flux densities is a more generic problem. Components B 
and C are clearly a merging pair of images associated with the tangential critical curve, and in 
most lens systems would be the brightest images. To reproduce the observed flux ordering in 
B1359-I-154, the critical curve must pass closer to A than either B or C, so that A will be the most 
magnified .P^ The 3SIS model is unable to accomplish this, as A sits quite far from the tangential 
critical curve (Fig. 5a). Accounting for the relative fluxes will therefore require models with more 
free parameters. 



4.2. More Realistic Models 

We now attempt to improve the fits to B1359-I-154 by investigating more realistic lens models. 
The mass distributions of lensing galaxies are ellipsoidal rather than spherical, so we consider 
models with 3 singular isothermal ellipsoids (SIE) in an external shear field. The SIE is described 



It is conceivable that the brightness of A may be produced by microlensing due to substructure in the lensing 
mass distribution (Schneider & Mao 1998). However, the large required amplification makes this less likely. 
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Fig. 5. — Best-fit 3SIS model, (a) Critical curves. Observed (model-predicted) image and galaxy 
positions arc marked as crosses (squares). Note that A is far from the tangential critical curve, and 
thus cannot be highly magnified in this model, (b) Caustic curves. Prominent imaging regions are 
marked as A/B, where A is the total number of images and B is the observable number of images. 
The source sits at the intersection of the dotted lines, near a swallowtail catastrophe, (c) Contours 
of constant time delay, in steps of 3h~^ days outward from A (thin lines), except for the critical 
contours (bold lines) passing through the saddle points in this surface (B, D, E and F). Images 
trapped in the galaxy cores correspond to local maxima. 
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by the scaled surface mass density 



K{x,y) 



6/2 



(1) 



[(l-e)a:2 + (l + e)y2]i/2 



where 6 is the critical radius and e = (1 — /^)/(l + f '^) is the ellipticity, with / as the mass axial 
ratio. In one set of trials we fix the orientations of G, G' and G" according to the de Vaucouleurs 
surface brightness profiles listed in Table 2. We refer to this as the 3SIE/FIX model. However, 
because these orientations are not very well constrained by HST, we also run trials in which the 
position angles of the galaxies arc left as free parameters (3SIE model). Finally, we recognize that 
galaxies forming the core of a compact group may be truncated on the scale of the inter-galaxy 
separation. We thus consider pseudo-Jaffe (PJ) models (Keeton Sz Kochanek 1998), which have a 
truncation radius at which the density breaks from p cx 1/r^ to 1/r^. The PJ profile is described 



We set a = (X'3, roughly half the inter-galaxy separation in this system, and investigate models 
with spherical (3PJS) and elliptical (3PJE/FIX, 3PJE) pseudo-Jaffe deflectors. 

The 3SIE/FIX model has a fit statistic of Xtot = ^16 for NDF = 7, significantly improving 
upon the 3SIS model with regard to the image positions {Xpos — '^^^ densities are still 
poorly matched to their observed values {Xfix — 182). The model-predicted flux densities of A-C 
are each much fainter than observed, and the fit to image A is particularly poor (predicted 2.8 
mJy, Xfix A ~ ''"'') • -'-t is unlikely that any fine-tuning of the model will remove this discrepancy so 
long as the galaxy position angles remain fixed. 

The 3SIE model relaxes constraints on the galaxy orientations and provides the best fit of 
any model we investigate, Xtot ~ ^"^ NDF = 4. The critical curves, caustics and time delay 
surface are plotted in Fig. 6. In this model the critical curve passes closer to A than either B or C 
(Fig. 6a), ensuring that A will be the brightest image. The key ingredient here is the placement 
of the unlensed source, which lies very close to the cusp of a swallowtail catastrophe inside the 
tangential caustic (Fig. 6b). As the source approaches this cusp, A diverges in magnification and 
would split into three adjacent images upon crossing into the 11/8 region. Therefore A can be 
highly magnified by placing a source just outside the swallowtail cusp. The largest contribution 
to the 3SIE fit statistic comes from the flux density of C (predicted 10.6 mJy, XfixC ~ ^^)- 
This is not a serious concern, as a slight shifting of the critical curve from B toward C could 
reproduce the proper flux ordering, and may be accomplished through unmodelable substructure 
in the lensing mass distribution. The second largest contribution is from the RA of galaxy G 
{XgaiG ~ 15.6), which accounts for nearly all of Xgai — 16.0. This is exacerbated by the small 
estimated uncertainty (0'.'008) of this coordinate. If we have underestimated the error, or had used 
circular error bars of radius equal to the major axis of the error ellipse (0'.'024 for G), x^al would 
be much smaller (~ 2). 



by 



K{x,y) = 



6/2 6/2 



(2) 



[(1 - e)x^ + (1 + e)y2]i/2 [(i _ e)x^ + (i + e)y2 + (j2]i/2 



-15- 




Fig. 6. — Best-fit 3SIE model. Plots are the same as in Fig. 5. Note that the critical curve (a) passes 
closer to A than either B or C, ensuring that A is the brightest image. The high magnification of 
A is due to the position of the unlensed source, which lies just outside a swallowtail cusp (b). 
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We now turn to the truncated mass distributions. The 3PJS model has a fit statistic of 

Xtot ~ ^3^5 ^ significant improvement over the 3SIS model but with the same number of free 
parameters (NDF = 4). The flux density of A is 11.2 mJy compared to the observed value of 22.6 
mJy, but overall the flux densities are fit rather well {Xfix = ^6). The fit statistic is dominated by 
offsets from the galaxy positions {x^al ~ 290). The 3PJE/FIX model also does better than its 
SIE counterpart, offering a fit of Xtot — ^'^^ NDF = 7. This is dominated by the fit to the flux 
density of A (predicted 8.8 mJy, Xfix A — '^^)- '^^^ 3PJE model offers a best-fit of Xtot — ^5, but 
again, A is not the brightest image (predicted 11.3 mJy, Xfix A ~ 2^)- interesting to note 

that while the 3SIE model more accurately reproduces the flux density of A, the predicted galaxy 
position angles in the 3PJE model (69°, 18°, 62° for G, G' and G", respectively) are much closer 
to those of the observed surface brightness distributions (48°, 29°, 63°, respectively). However, 
the external shears in the 3PJE and 3PJE/FIX models (7 2± 0.25) are significantly larger than 
those in the 3SIE and 3SIE/FIX models (7 ~ 0.15). 

To test whether the large external shears predicted by each of our lens models may be due 
to a dark matter halo associated with the galaxy group, we also investigated models with four 
distinct mass distributions. The external shear was set to zero in these trials. Although one might 
expect the three primary lenses to define the group center, we were unable to find reasonable 
models in which the group halo resides within the Einstein ring. Problems arise because a fourth 
mass distribution that is even moderately concentrated significantly disturbs the caustic structure 
of the system, unless it is centered on one of the lensing galaxies. This makes it very difficult to 
account for the image properties without creating additional lensed components. Not surprisingly, 
we find that in the vast majority of our trials the group halo is pushed out of the Einstein ring 
toward the east, where it becomes a large external shear contributor. 

We can check whether the mass ratios of the galaxies, as estimated from their lens model 
critical radii b, are consistent with the observed luminosity ratios. In standard lens models and 
observations of lenses, b oc L^/^ (see Keeton, Kochanek &; Falco 1998), so we predict critical radius 
ratios of 1.0:0.62 ± 0.2:0.63 ± 0.2 for G:G':G" based on the I-band luminosities in Table 2. The 
model ratios - 1.0:0.59:0.86 (3SIS), 1.0:0.62:0.75 (3SIE/FIX), 1.0:0.80:0.82 (3SIE), 1.0:0.85:0.55 
(3PJS), 1.0:0.68:0.60 (3PJE/FIX), and 1.0:0.57:0.68 (3PJE) - are remarkably consistent with 
these predictions. This provides additional evidence that all three lensing galaxies are important 
to the model and indicates that the B-band (for ~ 1) mass-to- light ratios of the three lenses 
must be very similar. The compatible galaxy colors and the consistency of our fundamental plane 
redshift estimates further reinforce this argument. 

5. Discussion 

HST and VLBA observations demonstrate that the gravitational lens B1359-I-154 consists 
of six images of a single background radio source. This is the first example of a galaxy-scale 
lens system with more than four images. The unique configuration is produced by the complex 
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mass distribution of the lens, a compact group of three ~ 1 galaxies lying on the vertices of a 
triangle separated by (/.'T ~ 4/i~^kpc, inside the l'/7 diameter Einstein ring defined by the radio 
components. The outer images (A-D) have the morphology of a standard four-image lens, but the 
triangle of galaxies produces additional extrema in the gravitational potential at which two central 
images (E and F) form. Simple lens models consisting of three deflectors constrained by HST 
galaxy coordinates and VLBA radio component data naturally produce six images at the observed 
positions, but require finely-tuned ellipticities to account for the flux density ratios. While it 
would be premature to claim that we have accurately represented the lensing mass distribution 
(xL/NDF > 10 for all models), our initial modeling efforts demonstrate that B1359+154 can be 
explained as a true six-image lens system produced by three lensing galaxies. The consistency 
between the model-predicted deflector mass ratios and the observed luminosity ratios of the lens 
galaxies adds credibility to our results. 

The three primary lens galaxies appear to constitute the core of a compact group. The most 
concentrated groups with small numbers of galaxies are the Hickson groups, a quarter of which 
have 3 members (Hickson et al. 1992). Separated by only ~ Ah~^ kpc, the three galaxies are close 

enough so that they might be expected to be consumed in a merger on a time scale of ~ 0.1 
(Barnes 1985). Thus the incidence of high multiplicity arcsecond-scale lens systems may allow for 
an estimate of the merger rate of massive galaxies at z ~ 0.5 — 1.5. 

The relative simplicity of B1359+154 compared to the lenses found in the cores of rich clusters 
should make this system an excellent tool for studying galaxy halos in dense environments, and 
the relationship between the galaxies and their parent halos. One issue worthy of investigation is 
whether the galaxy mass distributions are truncated on the scale of inter-galaxy separation. The 
improved fits of the 3PJS and 3PJE/FIX models relative to the 3SIS and 3SIE/FIX models offer 
tentative arguments to this effect. Because predicted time delays in the PJ models are significantly 
higher than those in the SIE models, measured delays may be used to discriminate between the 
profiles. A second interesting issue is the placement of the group halo relative to the lensing 
galaxies. The image positions and flux densities can be reproduced by three mass distributions 
associated with the observed galaxies, so there is little need to add a fourth independent mass 
distribution inside the Einstein ring based on modeling arguments. A fourth deflector above the 
critical density for multiple imaging would place an additional radial caustic onto the source 
plane and disturb the delicate balance needed to account for the data without creating additional 
images. Thus if this halo exists within the Einstein ring, it is likely to be either subcritical (large 
core radius) or centered on one of the galaxies. However, the similar model-predicted mass-to-light 
ratios of the galaxies argue against the latter possibility. The large external shear fields (7 ^ 0.15) 
required by all of our models suggests that a group halo may be displaced from the primary lenses. 
It is interesting to note that the orientation of the shear axis is remarkably constant, running 
nearly east to west regardless of the deflector model employed. 

The six lensed components allow for more complicated models than we have investigated 
here, but the development of a truly robust model of the lensing mass distribution requires the 
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acquisition of additional constraints. The detection of lensed extended emission from the quasar 

host galaxy in both the HST optical and CFHT infrared images is particularly promising. Deep, 
high-resolution HST imaging using either WFPC2 or NICMOS could investigate the properties 
of the Einstein ring in detail and provide vital new constraints on the potential structure (see 
Kochanek et al. 2000b). Furthermore, if subcomponents of D, E and F can be detected using deep 
global VLBI imaging, the relative orientations of the radio substructure could place additional 
constraints on the local magnification matrices. Finally, the measurement of differential time 
delays would also offer important constraints on the mass model, as these directly probe the 
lensing potential at the image positions. One might expect the source to be variable, based on 
the compactness of the radio components, though there is no evidence for significant variability 
at this point. Radio snapshots of B1359+154 should be routinely obtained to determine whether 
the variability is large enough to undertake a monitoring program using current instruments. 
Otherwise, the enhanced VLA may be able to determine time delays from the component light 
curves even if the source variability is only at the few percent level. 

It seems unlikely that the lens is a chance projection of three galaxies (an a posteriori 
probability of ~ 10~^ in a sample of 60 lenses), and their statistically identical colors further argues 
against this possibility. While the Next Generation Space telescope could directly determine the 
redshifts of all three galaxies, current telescopes can only measure the average redshift of the 
lenses unless they have narrow emission lines that are not apparent in the Myers et al. (1999) 
spectrum. Two possibilities are to search for Mg II metal absorption lines in the source spectrum 
or HI absorption features in the radio continuum once the average redshift is known. For — 1 
the Mg II absorption features would lie in the V-band where the source flux is 23 V mag, so the 
observation is difficult. Many groups have significant HI masses (Hickson 1997) and the radio 
flux density of ~ 90 mJy at 0.7 GHz (based on the 1.7 GHz VLBA flux density and a spectral 
index of a ~ 0.3 at low frequency) may be high enough to search for absorption features. X-ray 
observations are less promising because very deep observations would be required to detect the 
group (a bright 10^^ erg s~^ group at = 1 has an X-ray flux of only 4 x 10~^^ ergs cm~^ s~^) 
and the high resolution of the Chandra Observatory would be needed to distinguish emission from 
the group and the lensed sources. 
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A. The Statistics of Perturbing Galaxies 



Here we expand on the discussions of the statistics of perturbing galaxies of gravitational 
lenses in Keeton, Kochanek Sz Seljak (1997) and Kochanek & Apostolakis (1988) to include the 
effects of the galaxy luminosity function. We also emphasize the important distinction between 
a primary lens galaxy possessing a nearby neighbor and that neighbor significantly perturbing 
either the lens model or the caustic structure. We consider only galaxies which are clustered with 
the primary lens galaxy. Galaxies and clusters along the line of sight are an important source 
of weak shear perturbations but clustered galaxies dominate the perturbations (see Keeton et 
al. 1997; Barkana 1996). A key distinction to bear in mind is that all lens galaxies will have 
close neighbors (or satellites) if examined closely, so it is only the degree to which they modify 
the caustic structure or perturb the images of the lens which is relevant. Moreover, the 0.1 — 1% 
astrometric precision of most lens observations permits the detection of very small perturbations 
to lens models unless they are degenerate with other properties of the mass model. The sensitivity 
of models to small perturbations means that neighbors may be important to all lens models while 
remaining unimportant for most statistical properties of lens systems (image separations, image 
multiplicities, cross sections, etc.).|^ 

Following Keeton et al. (1997) and Kochanek & Apostolakis (1988) we model the distribution 
of neighbors using the correlation function, S,{r) = {r/rQ)~^ with — 5h~^ Mpc and x — l-^S 
(e.g. Peebles 1993), which when projected onto the lens plane becomes C2{R) — 4rQi?^~'^. We 
assume that the correlation function can be extrapolated to the small physical scales relevant 
to the calculations. We will add the effects of the distribution of galaxies in luminosity by 
including a luminosity function, dn/dL = n,,(L/L*)" exp(— L/L*) where ~ 0.01(/iMpc)^^ 
and Q ~ — 1 (e.g. Kochanek 1996). We use an SIS lens model to estimate the critical radius 
h = h^{L / LifY/"^ Dls / Dos and the Faber-Jackson relation to convert between luminosity and 
velocity dispersion, L/L* = [a/a^)^. Dql, Dqs and D^s are standard angular diameter distances 
and 5* = 47r((T*/c)^ = 1'.'45 for early-type lens galaxies (Kochanek et al. 2000a). The characteristic 
scale for finding neighboring galaxies is given by the probability of finding a galaxy with a 
comoving density n^^ inside the critical radius 6# of an lens galaxy. 



To = 27rn* / RdRi2{R) ^ 20n,r^{KDoLf-^ (Al) 

Jo 

hDoL 



0.005 



0.01(/iMpc)- 



n 1.75 



5/i-iMpc 



5h ^kpc 



1.25 



(A2) 



We will discuss only pairs of galaxies. While the probability of three galaxies can be determined from the 
three-point correlation function, it has never been measured on the kpc scales corresponding to the separations of 
galaxies in compact groups. An empirical determination based on a local redshift survey would probably be the most 
accurate approach to making the estimate. 
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The probability of finding a galaxy brighter than luminosity L is 



CXD 



t(> L) = To 



(iAexp(-A)A' 



a 



(A3) 



where \ = L/L^. The integral logarithmically diverges as L ^ for a = — 1, but the divergence is 
slow, as the integral is 0.2, 1.8 and 4.0 for L/L* = 1, 0.1 and 0.01, respectively. The result (A3) is 
an underestimate for several reasons. First, early- type galaxies are more strongly clustered than 
galaxies in general, so our value of rg is low. Second, for detecting neighbors we should include 
both early-type and late-type galaxies, and the total density of all types of galaxies is larger than 
0.01(/i-^Mpc)-3. 

The third problem is that the presence of multiple galaxies modifies the cross sections for 
finding the system as a gravitational lens. Kochanek &; Apostolakis (1988) showed that for galaxies 
with equal critical radii ha separated by distance I, the total cross-section is approximately 27r6o 
for t > 4bo, irbl [2 - (2 - l/2bo)f for 2bo < £ < 4bo and 7r(26o -^/2)2 for i < 2bo. The cross-section 
has a minimum at i = 2bo where the centers of the two lenses are deflected to the same point on 
the source plane. Where the total cross-section is smallest, the cross-section for more than two 
images and the average image magnification are largest. In short, as the cross-section is dropping, 
the amount of magnification bias will rise in compensation. For two galaxies with critical radii 6i 
and 62 separated by i we can derive several general results. There are no critical lines on the axis 
separating the two galaxies once £ < {b^^ + ^2^^)^' while the lens can generate images bracketing 
both lenses once i < 2{bi + 62)- If the separation is smaller than i < 2max(6i,&2) then it is 
possible to find solutions with no images in between the two lenses. Kochanek ct al. (2000b) used 
this to find a model for the four-image lens B1608+656 using two singular galaxies. The total 
multiple image cross-section is 7r(6f -|- 62) when the lenses are well separated, rising to 7r(6i -|- 62)^ 
when they are perfectly aligned. The dependence of the cross section on separation should be 
similar to the equal critical radius case, with a minimum total cross-section and a maximum large 
image multiplicity cross-section at an intermediate separation £ ~ (61 -|- 62). Because b oc L^/^, 
relatively low luminosity galaxies can produce substantial cross-section enhancements - a and 
O.IL^ galaxy have a cross-section when aligned which is 60% higher than that for the two galaxies 
separately. We can estimate the magnitude of the effect by computing the fractional "excess" 
cross-section 27r6i62/7r6i for aligned lenses over well separated lenses, assuming that it holds for all 
lenses separated by less than £ < 2{bi + 62)- Defining u = 61/6* the fractional excess cross-section 
from less massive neighboring galaxies (62 < &i) is of order 



For u = 1 {bi = 6*) the excess cross-section is IItq ~ 6%, while for u = 0.1 it is only O.Otq. 
While a numerical calculation including the effects of magnification bias is needed to obtain the 
exact impact of neighbors on lensing statistics, this order of magnitude calculation shows that the 
enhancements are modest. 




2 



(A4) 
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Keeton et al. (1997) considered a very simple model for the importance of neighboring galaxies 
by estimating the shear that a companion can produce at the center of the primary lens galaxy. A 
neighbor of luminosity L located at projected radius R produces a shear 7 = (h^DoL/'^R){L/ L^Y^"^ 
at the center of the primary lens. The probability of finding a shear larger than 7 is 

r(> 7) ~ 0.4ro7-^-^^ / dAA^+^^'^^Z^ exp(-A) (A5) 
Jo 

where A = L/L^,. The integral converges (for standard parameters q + (3 — x)/2 — —0.38 and the 
value of the integral is 1.4), and the perturbations are weakly dominated by low luminosity galaxies. 
While the probability of a lens having a neighbor diverges, the probability that the neighbor 
produces a large perturbation on the lens is finite and relatively small, r(> 7) ~ 0.6to/7^'^^. This 
estimate of the perturbations makes little sense for perturbing galaxies inside the Einstein ring or 
for estimating the perturbations to images at the Einstein radius. If wc instead estimate the shear 
produced at the closest point on the unperturbed critical line of radius 60 of the primary lens, we 
find 

r(> 7) ^ TO dAA" exp(-A) 

where the + (— ) branch corresponds to perturbations from outside (inside) the ring. The integral 
is not analytic, but under the approximation that |1 — (1 ± x)"| = we obtain the same result as 
in A3. The Keeton et al. (1997) approximation underestimates the total perturbation by about 

a factor of 2. The new approximation still overestimates the contribution from images inside the 
Einstein ring, because a perturbing galaxy located at the center is counted as producing a shear 
perturbation at the Einstein ring when its effects are indistinguishable from the monopole of the 
primary lens. 



1- 1± 



2760 



1.25 



(A6) 
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Table 1. VLB A Data 



Comp ARA (") ADec (") 5i.7 (mjy) b^aj (mas) af ^ 



Al 


= 


= 


22.6 


9 


-0.29 


Bl 


-0.490398 ± 0.00003 


-1.252428 ± 0.00003 


14.5 


6 


-0.32 


CI 


-0.311420 ±0.00002 


-1.669560 ± 0.00002 


19.8 


6 


-0.27 


Dl 


+0.962772 ± 0.00010 


-1.368612 ± 0.00010 


4.5 


5 


-0.32 


El 


+0.608880 ± 0.00008 


-1.142871 + 0.00008 


6.0 


3 


-0.23 


Fl 


+0.422273 ± 0.00025 


-0.963759 ± 0.00025 


1.8 


3 


-0.38 


A2 


-0.013783 ± 0.00013 


-0.006793 ± 0.00013 


3.3 


16 




B2 


-0.481571 ± 0.00032 


-1.292351 ± 0.00032 


1.4 


8 




C2 


-0.340726 ± 0.00034 


-1.622677 + 0.00034 


1.3 


7 




Note. 


— Model-fit components 


for B1359+154 1.7 GHz 


VLBA data (combined epochs). 


Positions 


are offset from RA 14** 


01™ 35^5476, Dec +15' 


13" 


25"646 (J2000). 


Nominal 



uncertainties on the flux densities are tiircc times tfie map rms of 45 /ijy/beam. Positional 
uncertainties are taken as the beam size (~ 10 mas) divided by the signal-to-noise for each 
component. The major axes of the gaussian components are given by bmaj- Spectral indices 
J between 1.7 and 5 GHz arc calculated using the combined flux density of the primary and 
subcomponents, where applicable. The 5 GHz data is from Rusin et al. (2000a). Assuming 5% 
errors on both the 1.7 and 5 GHz flux densities due to vajriability and model-fitting, the error 
on the spectral indices is ~ 0.07. 
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Table 2. HST Astrometry and Photometry 



ID 


ARA (") 


ADec (") 


I (mag) 


V-I (mag) 


Re (") 


e 


PA (°) 


A 


=0 


=0 


24.01 ±0.06 


0.68 ±0.12 








B 


-0.483±0.007 


-1.253±0.009 


24.29 ±0.11 


0.43 ±0.15 








C 


-0.323±0.007 


-1.640±0.003 


24.01 ±0.05 


0.27 ±0.13 








D 


0.957±0.008 


-1.357±0.008 


25.19 ±0.08 


0.61 ±0.17 








E 


0.627±0.013 


-1.129±0.011 


25.12 ±0.17 


0.45 ±0.19 








F 


0.426±0.017 


-0.951±0.028 


26.25 ± 0.23 


0.29 ±0.23 








G 


0.160±0.008 


-0.935±0.024 


22.68 ±0.28 


3.31 ±0.55 


0.71±0.29 


0.11±0.23 


48±55 


G' 


0.658±0.005 


0.690±0.008 


23.69 ± 0.24 


3.13 ±0.29 


0.16±0.06 


0.69±0.14 


29±8 


G" 


0.r)75=().0l)7 


J..r,l(i±O.I)ll 


2:170 = ().:-!:-! 


:-!.2:-l = 1.01 


0.1 1±0.()7 


1). 19 = 1). 20 


(i:-!±7 



Note. — The three lensing galaxies arc fit to dc Vaucoulcurs profiles. The effective radii {Re), ellipticity 
parameters e = 1 — / and position angles (PA) of these fits are listed. 
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Table 3. Nearby Objects 



ID 


ARA (") 


ADer ("\ 


I 


V-I 


I-K 


1 


6j 


Notes 


G 


0.2 


-0.9 


22.68 ± 0.55 


3.31 ±0.97 


~2.0 






K3, Group 


G' 


0.7 


-0.7 


23.69 ±0.14 


3.13 ±0.35 


~3.0 






Kl, Group 


G" 


0.6 


-1.5 


23.70 ±0.12 


3.23 ± 1.03 


~3.0 






K2, Group 


SI 


4.4 


9.3 


24.39 ± 0.11 


> 1.74 










S2 


11.5 


-2.9 


24.94 ±0.12 


> 1.43 










Gl 


-13.5 


-7.5 


19.08 ±0.10 


1.36 ± 0.03 




0.076 


64 




G2 


-16.8 


1.7 


22.73 ±0.10 


1.42 ± 0.05 




0.013 


-81 




G3 


0.6 


-10.6 


23.30 ±0.10 


0.79 ± 0.04 




0.017 


-3 




G4 


2.1 


5.2 


23.52 ±0.11 


0.94 ± 0.06 




0.023 


18 




G5 


3.7 


-1.3 


23.52 ±0.11 


> 2.41 


~4.4 


0.042 


-85 


K4, Group 


G6 


-17.5 


6.7 


23.58 ±0.11 


1.09 ± 0.06 




0.008 


-67 




G7 


-9.1 


-5.1 


23.97 ±0.11 


1.17 ±0.06 




0.012 


66 




G8 


-7.1 


8.4 


24.31 ±0.11 


1.48 ± 0.08 




0.009 


-38 




G9 


-1.8 


1.5 


24.50 ± 0.12 


> 1.37 


~3.5 


0.030 


-38 


K5, Group 


GIO 


6.8 


-12.2 


24.70 ±0.12 


0.63 ± 0.09 




0.007 


-31 




Gil 


12.2 


-2.6 


24.80 ±0.12 


1.54 ± 0.09 




0.007 


-82 




K6 


-6.5 


-6.9 


> 25.0 




>4 






K6 


Tl 


-32.3 


-67.8 


18.26 ± 0.10 


1.38 ± 0.03 




0.023 


26 




T2 


-47.0 


-9.9 


18.85 ±0.10 


2.20 ± 0.03 




0.027 


79 




T3 


-34.0 


-50.1 


19.23 ±0.10 


2.12 ±0.03 




0.018 


35 





Note. — SExtractor position offsets from A are given in arcscconds. G,G',G" photometry and positions are from 
component fits in Table 2. Galaxies on the WFPC2 image within 20" of G axe labeled Gx, and those detected only 
in the Rusin et al. (2000a) infrared observations are labeled Ka;. Stars are labeled as Sx. Galaxies outside of 20" 
with large estimated shears are included as Tx. Tidal shear estimates 7 assume SIS halos with the same redshift and 
mass-to- light ratio as G, and the shear position angles 6~i are relative to G in degrees counter-clockwise from North 
(see Lehar et al. 2000). 
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Table 4. Lens Models 



Type 
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b 


Xl 


yi 


e,7 


PA 


Img 


AtiA 












(") 


(") 


(") 




n 




(/i~^days) 








QQTG 




n Ar\r\ 
U.4UU 


n 1 70 
u.i /y 


— u.yuo 


n nnn 
U.UUu 


r\ n 
U.U 


r> 


Q 7 

y. / 


Xtot 


-1 CfJQ 


Q 

.O 




Vjr 


u.zou 


u.uou 




nnn 

U.UUU 


U.U 


n 

\^ 


y .o 




- 

Apnt 


1 ^41 fi 




vjr 




U.OO ( 


— i.ooy 


n nnn 

U.UUU 


n n 

U.U 


u 


OO.O 




\flx 


154.6 




olltJcLl 








U. lO^ 


87 Q 
o / .y 


Hi 


OU.O 




Xgal 


73.6 




src 




u.oou 










ou.y 








QQTT? /TT'TV 


Vjr 


n /I no 


u.ioo 


— u.yiz 


n 1 7n 
u.i ( u 


/IS n 


"D 


Q fi 

y.o 


^2 
A.tot 


= 215.i 


i 










—0 71 S 

L». 1 -LO 
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fi^ n 
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T? 
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04.0 
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16.0 
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1 7 


^2 

A-tot 


= 537.9 




G' 


1.170 


0.618 


-0.693 


0.000 


0.0 


c 


16.6 




Xpnt 


202.1 




G" 


0.754 


0.649 


-1.416 


0.000 


0.0 


D 


66.5 




X%x = 


45.8 




shear 








0.235 


96.2 


E 


69.3 




Xgal 


290.0 




src 




0.223 


-0.919 






F 


71.2 








3PJE/FIX 


G 


1.348 


0.109 


-0.964 


0.171 


48.0 


B 


20.5 


Xtoi 


= 99.8 






G' 


0.914 


0.664 


-0.700 


0.448 


29.0 


C 


19.7 




Xpnt 


3.9 




G" 


0.812 


0.581 


-1.488 


0.506 


63.0 


D 


66.1 




Xflx 


44.0 




shear 








0.255 


94.0 


E 


66.7 




Xgal 


51.9 




src 




0.242 


-0.941 






F 


68.2 








3PJE 


G 


1.095 


0.134 


-0.962 


0.581 


69.2 


B 


27.5 


Xtot 


= 64.7 






G' 


0.626 


0.658 


-0.693 


0.665 


18.4 


C 


25.8 




Xpnt 


1.1 




G" 


0.745 


0.580 


-1.489 


0.603 


62.0 


D 


73.9 




Xflx — 


45.1 




shear 








0.262 


95.8 


E 


73.6 




Xgal 


18.5 




src 




0.230 


-0.929 






F 


74.9 









Note. — The five lines for each model give the properties of the lensing galaxies (G, G' and G"), external shear, 
and recovered source coordinates: 6 is the critical radius, xi and yi are the lens coordinates, e = 1 — / where / is 
the axis ratio of the ellipsoid, 7 is the strength of the external shear, and PA is the major axis position angle of the 
ellipsoid or the direction (to a tidal perturbation) producing the shear. The time delays for each imago arc defined 
relative to image A as AtiA = U—tA and were computed for an f2o = 0.3 fiat universe and a lens redshift Zi = 1. The 
final column gives Xto* and the contributions from the astrometry of the images (xpos), the fiux ratios of the images 
{Xfix) and the positions of the lens galaxies (Xgai)- The models assume that the uncertainties in the positions of the 
radio sources are 1 mas and that the fractional uncertainties in their fiux densities are 10%. 



